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Waves propagating in the relativistic electron-positron or ions plasma are investigated in a frame 
of two-fluid equations using the 3 + 1 formalism of general relativity developed by Thorne, Price 
and Macdonald (TPM). The plasma is assumed to be freefalling in the radial direction toward 
the event horizon due to the strong gravitational field of a Schwarzschild black hole. The local 
dispersion relations for transverse and longitudinal waves have been derived, in analogy with the 
special relativistic formulation as explained in an earlier paper, to take account of relativistic effects 
due to the event horizon using WKB approximation. PACS: 95.30.Qd, 95.30.Sf, 97.60.Lf 
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1. Introduction 

In recent year plasma equations and general relativ- 
ity are usually considered together. The Coulomb 
potential of charge particles due to coupling is much 
stronger than the gravitational potential and are of- 
ten neglected in the Newtonian approximation. But 
the mean gravitational field for certain astronomi- 
cal objects like galactic nuclei or black holes may 
be strong and the observation of magnetic fields in- 
dicates that a combination of general relativity and 
plasma physics at least on the level of a fluid descrip- 
tion is appropriate. The plasma in the black hole en- 
vironment may act as a fluid and black holes greatly 
affect the surrounding plasma medium (which is 
highly magnetized) with their enormous gravitational 
fields. Hence plasma physics in the vicinity of a 
black hole has become a subject of great interest in 
astrophysics. It is therefore of interest to formulate 
plasma physics problems in the context of general 
relativity. 

Thorne and MacDonald [TJ [5] have introduced 
Maxwell's equations in 3+1 coordinates, which pro- 
vides a foundation for formulation of a general rel- 
ativistic (GR) set of plasma physics equations in 
the strong gravitational field of both the nonrotat- 
ing and rotating black holes and the "membrane 
paradigm" [3] is a good example of such a formal- 
ism in which the electromagnetic equations and the 
plasma physics at least look somewhat similar to the 
usual formulations in flat spacetime while taking ac- 
curate account of general relativistic effects such as 
curvature. The membrane paradigm is mathemat- 
ically equivalent to the standard, full general rela- 
tivistic theory of black holes, so far as all physics 
outside the horizon is taken into account. 

Sakai and Kawata (SK) [4. have developed the 
linearized treatment of plasma waves using special 
relativistic formulation. Such an investigation of 
wave propagation in a general relativistic two-fluid 
plasmas near a black hole is important for an un- 
derstanding of plasma processes. That is, what hap- 



pens when the plasma are assumed to be freefalling 
onto the black hole. The study of plasma wave in 
the presence of strong gravitational fields using the 
3 + 1 approach is still in its early stages. Zhang [5j[6] 
has considered the care of ideal magneto hydrody- 
namics waves near a Kerr black hole, accreting for 
the effects of the holes angular momentum but ig- 
noring the effects due to the black hole horizon. Hol- 
comb and Tajima [7], Holcomb [8], and Dettmann 
et. al. 9 have considered some properties of wave 
propagation in a Friedmann universe. Daniel and 
Tajima [10] studied the physics of high frequency 
electromagnetic waves in a strong Schwarzschild plasma. 
Marklund et al [H] have found a mode representing 
high frequency plasma oscillation in a charged two- 
component plasma using the exact 1+3 covariant 
dynamical fluid equations in the presence of elec- 
tromagnetic fields about a Friedmann-Robertson- 
Walker model by ignoring the fluid's thermal ef- 
fects. Servin et al [12] and Kleidis et al [13] H4] 
have studied the propagation of gravitational waves 
in a collisionless plasma with an external magnetic 
field parallel to the direction of propagation, while 
Forsberg et al [15] have presented an investigation 
of nonlinear interactions between gravitational ra- 
diation and modified Alfven modes in astrophysical 
dusty plasmas. 

There is also work on fluid dynamics and kinetic 
gas theory in the context of cosmology. The book by 
Bernstein [16] treats gas kinetics in the Friedmann- 
Lemaitre-Robertson- Walker (FLRW) model. How- 
ever, there are relatively few relativistic cosmologi- 
cal investigations that take into account plasma ef- 
fects and the behavior of matter in the presence of 
electromagnetic fields [HI [HI [El EDJ EH E21 E3]. 
Therefore, the general relativistic treatment of plas- 
mas, both in astrophysics as well as in plasma physics, 
seems to be a field open to investigation. 

A plasma can propagate both linear and non- 
linear waves. Linear refers to the simplifying ap- 
proximations that are possible for small amplitude 



waves like, Alven and high frequency electromag- 
netic waves, and nonlinear refers to large amplitude 
phenomena not predicted by linear models. In this 
paper, the set of two-fluid equations for collision- 
less ideal plasma are used to make an initial at- 
tempt to be the discovery of an instability caused 
by the general relativistic term in the dispersion 
relations for transverse (electromagnetic) and lon- 
gitudinal (electrostatic) waves using action princi- 
ple for a hot plasma developed by Hcintzmann and 
Novello [24]. Similar multifluid equations have re- 
cently been used to calculate local dispersion laws 
for plasma waves in strong and weak gravitational 
fields; see, e.g., Buzzi et al [251 EH], and Rahman et 
al [2H HH1 12E IHD] ■ Thus, the present work will form 
the essential basis of nonlinear, more complicated, 
investigations. 

In the present paper Sec. 2 summarize the 3+1 
formulation of general relativity. In Sec. 3 we re- 
view the two-fluid plasmas governing equations in 
Schwarzschild coordinates. The transverse and lon- 
gitudinal parts are separated by introducing a new 
complex transverse fields and velocities using Rindler 
coordinates in Sec. 4. In Sec. 5 the two fluid equa- 
tions are linearized for wave propagation by giving 
a small perturbation to fields and fluid parameters. 
We discuss the way in which the unperturbed fields 
and fluid parameters and their derivatives with re- 
spect to z depend on the surface gravity of the black 
hole and freefall velocity in Sec. 6. In Sec. 7 the 
two-fluid equations are simplified using action prin- 
ciple developed by Heintzmann and Novello [23] . 
The dispersion relations for transverse and longitu- 
dinal waves are developed and solved using the an- 
alytical method developed by Mikhailovskii [31] in 
Sec. 8 and 9. Finally, we present our remarks in Sec. 
10. Here, we use units in which G = c = fcs = 1. 

2. 3+1 Formalism of Schwarzschild 
Spacetime 

Our work presented in this paper is based on the 
3+1 formulation of general relativity developed by 
Thome, Price, and Macdonald (TPM) [HUE]. The 
basic concept behind the 3 + 1 formulation of gen- 
eral relativity is to select a preferred set of space- 
like hypersurfaces which form the level surfaces of 
a congruence of timelike curves. A particular set 
of these hypersurfaces constitutes a time slicing of 
spacetime. The hypersurfaces considered here are of 
constant universal time t. In the 3 + 1 formulation, 
the Schwarzschild metric is given by 

ds 2 = -a 2 dt 2 + \dr 2 + r 2 (d9 2 + sin 2 ^ 2 ), (1) 



where the lapse function a is given by 



= y/l - 2M/r. 



(2) 



The hypersurfaces of constant universal time t de- 
fine an absolute three-dimensional space described 



by the metric 



1 



ds 2 = —dr 2 + r 2 (d0 2 + sin 2 6»^ 2 ). (3) 
or 

We consider a set of fiducial observers (FIDOs) 
having properties as mentioned by Buzzi et al [25] 
remaining at rest with respect to this absolute space. 
For a detailed concept of a set of fiducial observer 
(FIDOs) see the membrane paradigm book [3]. To 
make the local measurements of all physical quanti- 
ties FIDOs use a local Cartesian coordinate system 
that have basis vectors of unit length tangent to the 
coordinate lines: 
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(4) 



For a spacetime viewpoint rather than a 3 + 1 split 
of spacetime, the set of orthonormal vectors also 
includes the basis vector for the time coordinate is 
given by 

d - 1 9 fK\ 

66 " dr^ adt' [b) 
The lapse function a plays the role of a gravitational 
potential and thereby governs the ticking rates of 
clocks and redshifts. We can calculate the gravita- 
tional acceleration felt by a FIDO from the lapse 
function as follows jTJ [2j [3] : 



-Vina = o e fi 



1 M 

a r 2 



(6) 



Equation (jSJ) shows that far from the black hole 
event horizon the gravitational acceleration become 
weak and approaches the Newtonian value for flat 
spacetime. However, near the horizon, the gravi- 
tational acceleration approaches infinity as a — > 0. 
The rate of change of any scalar physical quantity or 
any three-dimensional vector or tensor, as measured 
by a FIDO, is defined by the derivative 



D (\ d 

= h v • V 

Dr \adt 



(7) 



v being the velocity of a fluid as measured locally 
by a FIDO. Since all the quantities are measured 
locally by the FIDO, all the vector quantities are 
neither covariant nor contravariant. 

3. Two-fluid Equations 

We consider two-component plasma such as an electron- 
positron or electron-ion. The basic equations (the 
equation of continuity and Maxwell's equations) for 
each fluid of species s with velocity v s , number den- 
sity n s , and the relativistic Lorentz factor 7 S as de- 
rived by TPM [Till [3] and Buzzi et al [25] [25] using 
the 3+1 formulation are given by 

d 

— (j s n s ) + V • (afy s n s v s ) = 0, (8) 
at 



and 



V B 

V • E 

dB 

~dt 
<9E 

~dt 



0. 

47TCT, 

-V x (aE), 

V x (qB) - 47raJ, 



where the charge and current densities are defined 
with charge q s of each fluid species by 



(13) 



where s is 1 for electrons and 2 for positrons (or 
ions). The presence of lapse function a signifies 
the general relativistic effect around a Schwarzschild 
black hole. 

Using Maxwell's Eqs. (PJ and ((EJ the equa- 
tions for the conservation of energy and momentum, 
as derived by TPM [TJHE] and Buzzi et al [25l [26] . 
for Schwarzschild black hole are given by 

= -V-S s + 2a-S s , (14) 

a at 

Id 1 <— > 

- — S s = e s a -V- aW s ), (15) 

a at a 

where the energy density e s , the momentum den- 
sity S g , and stress-energy tensor Wjj k for the elec- 
tromagnetic field as 



±(E 2 +B 2 ) 



S s = — E x B, 

s 4tt 



Wj> k = — (E 2 + B 2 W fe - — (E ] E k + BW k ). (16) 

8"7T 47T 

For a perfect relativistic fluid of species s in three- 
dimensions, the energy density e s , the momentum 
density S s , and stress-energy tensor W| fc correspond- 
ing to the above equations for electromagnetic field 
are 

e s = 7s 2 (e, + P s v 2 ), S s = 7 2 ( £s + P s )v s , 
^ = 7 2 ( £s +P s K:« s fc +P^ fc , (17) 

where v s is the fluid velocity, n s is the number den- 
sity, P s is the pressure, and e s is the total energy 
density defined by 



m s n s + P s /(7 9 - 1). 



(18) 



The gas constant 7 g take the value 4/3 for T — > 
oo and 5/3 for T — > 0. Using the conservation of 
entropy the equation of state can be expressed by 



D 



0. 



(19) 



The energy and momentum conservation Eqs. (|14l) 
and (fH)]) coupling with each single perfect fluid of 



species s to the electromagnetic field are 

^( es + p s )]_v.[ 7s 2 ( £ 



(9) 
(10) 


adt s adt 
+j s q s n s ~E ■ v s -( 


(11) 


^ s (e s + P s )(- 
\a 


(12) 




+v s x B) + v s 



+ P,)a- 
■ V 



7 s g s n s E ■ v s 



v s = 0, (20) 
f- VP S - ~f s q s n s (E 
1 9 



a dt 



P 



+7 2 (e s + P s )[v s (v s -a)-a]=0 



(21) 



Here, e s is the internal energy density and P s is the 
fluid pressure. Similar energy and momentum con- 
servation equations have previously been obtained 
by Buzzi et al [25l [26]. If, now, one sets a = 1 
so that the acceleration goes to zero, these equa- 
tions reduce to the corresponding special relativis- 
tic equations as given by SK [J]. The curl and di- 
vergence operators in the 3+1 set of equations are 
covariant and can be derived in locally cartesian co- 
ordinates instead of the spherical three-metric us- 
ing Rindler coordinate system, in which space is lo- 
cally Cartesian, provides a good approximation to 
the Schwarzschild metric near the event horizon in 
the form 



ds 2 



-a 2 dt 2 + dx 2 + dy 2 



dz 2 



(22) 



where 



2M, 



x = r H {6- w/2),y = r H 4>, z = 2r H (l - — )~ 1/2 . 

r 

. (23) 

The standard lapse function in Rindler coordinates 
becomes a = z/2ru, where th is the Schwarzschild 
radius. One of the advantages of the Rindler geome- 
try is that it gives an example of the essential ideas 
of the horizon and the 3+1 split without the dis- 
tracting complication of curved spatial three-metric. 
The transformation from the FIDO comoving (fluid) 
frame to a coordinate frame of the metric given in 
Eq. ([]} involves a boost velocity, which is a simple 
Lorenz boost with velocity Va in the radial direction 
defined by 

v ff =(l-a 2 )5. (24) 
Then the relativistic Lorentz factor becomes Tboost E 



(1 



-1/2 



I /a. 



4. Wave Propagation in Radial Di- 
rection 

Here we consider an incoming gravitational wave 
propagating in radial z direction toward the event 
horizon in presence of an external static magnetic 
field B = Po e i an d study gravitational waves exci- 
tation of small amplitude plasma waves, restricting 
our attention to a one-dimensional case. Introduc- 
ing the following complex variables, 



v s (z,t) = v sx (z,t) + iv sy (z,t), 
B{z,t) = B x (z,t)+iB y (z,t), 
E(z,t) = E x (z,t) + iE y (z,t), 



(25) 



the equation of continuity, Eq. ([8]), and Poisson's 
Eq. ([T0|) for each fluid species can be written as 



d d 

— (7 s n s ) + — (a^ s n s u s ) = 0, 
at dz 



dE z 
dz 



47r(gini7i + g 2 "-272)- 



(26) 
(27) 



Two transverse parts obtaining from adding the eg 
component multiplied by i to the e% component of 
the Maxwell's Eqs. jTTJ and ||ISJ| may be written in 
single form as 

l dz 2 dz dz dt 2 \dz J J 
d 

= 47rea— (n 2 7 2 v 2 - "i7i u ij- ( 28 ) 

The transverse and longitudinal parts of the mo- 
mentum conservation Eq. (|21|) can be separated 
out of the form 

Du s ( i A c>p s 



Dt 



dz 



i ap s 



+ (1 - t£)p»a - u s ( q s n sls -E . Vs + ) , (29) 



Dt 



q s n s j s (E - W S B Z + iu s B) 



( 1 dP s \ , s 

-u s v s p s a - v s yq s n s j s E ■ v s + — I , (30) 

where the suffix star is the complex conjugate, u s 
the z component of velocity, and the total energy 
density defined by p s = jg(e s + P„) — "f 2 (m s n s + 
T g P s ) with T g = 7 s /(7 9 — 1). In order to investi- 
gate the transverse electromagnetic waves it is more 
convenient to work from a combination of Eq. (|28p 
and momentum conservation Eq. (|29[) . The longi- 
tudinal waves can be investigated by combining the 
longitudinal components of the continuity Eq. (|26[) , 
Poisson Eq. (1271) . and the conservation of momen- 
tum Eq. (|50]). 

5. Linearized Equations 

We linearize the two-fluid equations by considering 
a small perturbation. We introduce the quantities 

n s (z,t) = rio s (z) + 5n s (z,t)), v s (z,t) = 8v„(z,t), 

P s {z,t) = P Qs (z)+SP s (z,t),-E(z,t) = SE{z,t), 

u s (z,t) = uq s (z) + 5u s (z,t), B(z,t) = 5B(z,t), 

Ps(z,t) = po s (z) + 6p s (z,t), 

B z (z,t) = B (z) + SB z (z,t). (31) 



The transverse set of two-fluid equations (|28|) and 
(|2T)1) are linearized using Eq. (|3ip of the form 



9 d 2 „ da d d 2 
1 3a 



dz 2 



dz dz 




= 47rea lno27o2-^ «oi7oi-^- ) , (32) 



d d 



auQ S — + — - u 0s a 
dz dt dz 



da iaq s j 0s n 0s B \ dSv, 



POs 



dt 



aq s jo s n 0s 

POs 



d dot 
au 0s — + — + u Qs — ) SE = 0. (33) 
dz dt dz ' 



The longitudinal set of equations (f2l)|) . (|2"7|) . and 
(|3"0|) are linearized to obtain 



7o s 

d 
dz 



d 



da 



du 



— + u 0s a— 
dt dz 

J {nQs"fosU 0s 



O.s 



Sn, 



UOs dt +a dz 



da 
dz 

dSE z 



1 dn 0s , _ 2 du o 
+ 37 0s ii 



n 0s dz 



5u s = 0, (34) 



(9- 



47re(n 02 7o 2 - U01701) + 47re(7 02 (5n 2 



-7oi(5ni) + 47re(n 02 Mo27o2^2 - UoiUoiToi^i)) (35) 



and 



2 / -i 2 \ 

77: + w o s o;— + 7 0s "( 1 + u 0s) 



d_ 

di 



aq s na s 

P0s70s 



5E Z 



u 0s a- 



d_ 

dz 
du 0s 



dz 



dz 

a dPo s 1 da 
p 0s dz j 2 s dz 



— 2 1 U 0s-^7 

7o s "Os I Ai* V at 



7o s 7g-Pos 

POs 



A). 

TgPos ( 1 dP 0s 

POs 



wo s 7o s a 



d 

+ a "H~ 

dz 

1 dn 0s 




0. (36) 



6. Dependence of the Unperturbed 
Values on z 

Since the plasma is assumed to be falling in radial 
direction, the infall velocity can be defined as 



u 0s (z) = Uff(z) = [1 - a 2 {z)}i. 



(37) 



The unperturbed number density, pressure, temper- 
ature, and magnetic field can be determined directly 
from the equation of continuity. From Eq. (1261) it 



follows that r aj 0s n 0s UQ S = const. = rj I aH"fH , nH u H, 
where the values with a subscript H are the limiting 
values at the event horizon. The freefall velocity at 
the horizon becomes unity so that uh = 1. Since 
uos = Vf£, 70s = and hence aj 0s = a H lH = 1. 
Also, because vg — (rjj/r)3, the number density, 
unperturbed pressure, temperature profile, and un- 
perturbed magnetic field for each species can be 
written as follows: 

n 0s (z) = n Hs v%(z), Pqs(z) = P ffs % 7s (z), 



Tqs — Pffs% 



3(7 9 -l) 



(z),), B (z) = B H v s (z), (38) 



with ks — 1 and Pq s = fcsno s To s . The derivatives 
of the above quantities with respect to z in Rindler 



coordinates expressed as 



dvg 
dz 
dn 0s 
dz 



du 



a 1 dBo 



dz 2th Vff ' dz 

3a n 0s dP 0s _ 
2ru Vg ' dz 



4a B 

3a "{ g Po s 
2r H vl 



(39) 



7. The WKB Approximation 

We consider the infinitesimal displacements of waves 
with small amplitude toward the horizon in WKB 
approximation. We can write all the perturbations 
of the form fo(z) exp(« J k(z)dz — iut), where fo{z) 
and k(z) are slowly varying function of z. Since the 
freefalling in any gravitational field occur only for lo- 
cally, the correct WKB solution can be derived from 
the action principle developed by Heintzmann and 
Novello [24]. From this standpoint, the local disper- 
sion relation and the instability can be anticipated 
from basic principle. The only scale in this problem 
is the black hole radius rjj, and the amplitude is 
small enough, we can ignore the internally reflated 
wave as long as rjj / A ^> 1 because the amplitude of 
this wave vanishes as e~ rff / A , where A = 2w/k(z). 
Thus, the set of transverse two-fluid equations, Eqs. 
(|32|) and ([33]) become 

(a 2 k 2 - cu 2 ) SE 

= i47reaw(no27o2^2 - "oiToi^i), (40) 

/ , . aq s j Qs n 0s B \ 

u) akuo s — oj H dv s 

V POs J 

-ia ™" 0s (aku 0s -u>)8E = 0. (41) 

POs 

Using Eq. (13T)1) . the set of longitudinal two-fluid 
Eqs. (j3"4l - (l3T>l) when Fourier transformed become 

n 0s 7o s (ak - u Qs uj) 6u s + (au 0s k - oj) Sn s = 0, (42) 



It is clear from the dispersion relations given in 
Eq. (|4"5]) that the general relativistic effects en- 
ter only in the ratio ui/a. If one uses local time 
dr = adt of the FIDO instead of a global coordi- 
nate time t, then these dispersion relations reduced 
to the special relativistic version as it should be ac- 
cording to the Einstein relativity principle. 

If we are looking the background of the infall 
radial velocity of the fluid species toward the event 
horizon which is near unity closed to the event hori- 
zon (i.e., at a — > ) and decreases with the distance 
from it to zero (i.e., at a — > 1 ) then the Eq. ([45]) has 
simple analytic solutions at ^ ^ kuQ Sl ^ » kuo s , 
and — <C kuos- Let us consider these solutions. 

8.1 Electromagnetic waves with — ^> kuo s 
We considered the first case corresponds to the waves 
having phase velocity larger then the infall radial ve- 
locity (i.e., at -*x 3> uq s ). In this case, the index of 
refraction can be approximated from Eq. (|45j) of the 
form 



2 2 
= 1-^ 



pi 



u p2 



1 



1 



(46) 



This equation shows that this waves have a reso- 
nence not for electron, depends on the combine os- 
cillations of electron with positron or ion. For electron- 
ion plasma, there is an asymmetry between the par- 
ticle species due to the small mass ratio m e /mi, giv- 
ing different order of magnitudes for the two local 
cyclotron frequencies, u) p i and uj P 2- Thus the reso- 
nances will occur vary different wave frequencies. 

For electron-positron plasma (m e = rrij, two 
species) Eq. (j45]) can be written with u) c \ = cu C 2 
and w 



pi 



^p2 



N 2 



(47) 



(auosk — bS) 5u s 



u Ts 



(ak — uq s uj) Sn s 



iaq s n 0s 



5E Z = 0, (43) 

POsJOs 

ikSE z = 47re(n 2702 - noiToi) + 47re(7 02 <5n.2 
-7oi<5ni) + 47re(n 02 M 27o2^ u 2 - n itioi7oi^ u i)- ( 44 ) 

8. Transverse Electromagnetic Os- 
cillations 

The dispersion relation for transverse electromag- 
netic waves may be obtained from Eqs. (|38[) and 
(1591 as 



J>_ _ UJ 2 pl (| - U 01 k) | 0J 2 2 (g - U 02 k) 

a 2 UQik-^-ui c i u a2 k-%+uj C 2 



(45) 



for either the electron-positron or electron-ion plasma. 
The local plasma frequency u ps = \J Ane 2 ~f' 2 s n 2 )s / p$ s 
depends on the local number density of each fluid 
species and the local cyclotron frequency 

co cs = e 2 jo s no s Bo/ pq s depends upon both the 
local number density and magnetic field. 



Which implies that two electromagnetic wave modes 
are exist: the upper branch represents the high fre- 
quency electromagnetic wave with frequency — > 
(2w 2 s + Lo 2 ^ 1 / 2 and another is the low frequency 
Alfven wave with frequency ^ < uj cs . On the other 
hand there appears a cut-off frequency ^ = (2uj 2 s + 
u^) 1 / 2 in the range of electromagnetic wave. A res- 
onance occurs when — = uo ps . In this approximation 
our results agree with the results found by Daniel 
and Tajima [10], and SK [4] corresponding to gen- 
eral relativity and ultrarelativistic limit. 

8.2 Electromagnetic Waves with — « kuo s 
When the phase velocity of the waves approaching 
the fluid's infall velocity, i.e., -^t ~ uq Si we obtain 
from the dispersion relation given in Eq. (|45|) , 



N 2 = 



ak 



(48) 



which shows that the electromagnetic waves reap- 
pear. Also, in this case we have Uq s = 1, i.e., the 



infall radial velocity is maximum and is known as 
the freefall velocity of the fluid species. 

8.3 Electromagnetic Waves with — <C kuo s 
Another solution can be simplified when the wave 
speed is far, or at least, is not very closed to the 
infall velocity (i.e., at <C uos) of the form 

(ttoifc - uj c i)(u 02 k - uj c2 ) = -Moiuoi^p! + w 2 2 ), 

(49) 

which for electron-positron plasma becomes 



x 0s 



(50) 



Since ^ -C kuo s < uj cs , electromagnetic waves in a 
strong magnetized plasma reappear at frequency be- 
low the cyclotron frequency. When the plasma den- 



sity is large and the field strength is small (2uj 2 s Uq s > 
uj 2 s ) the wave number becomes imaginary (fc 2 < 0) 
and the branch of a periodically damped and grow- 
ing (i.e., unstable) oscillations exist. Similar in- 
stabilities of a relativistic plasma were found by 
Mikhailovskii [31], and Zaslavskii and Moiseev [32] . 
Here the damping corresponds to Im(fc) > and 
growth to Im(fc) < 0. 

9. Longitudinal (Electrostatic) Os- 
cillations 

The longitudinal waves dispersion relation can be 
obtained from Eqs. ((42]), (g3]), and ((44]) of the form 



Kifc-^) 2 -^(fc-^oi) 2 



u p2 



(k 



(51) 



where v 2 -. 



^IglosPos/ Pos is the fluid's thermal 
velocity. Since the 7q s factor involved in the energy 
density po s in the denominator, the 7g s factor in the 
numerator cancels out and therefore, the thermal 
velocity of the two-fluid plasmas are frame indepen- 
dent. 

The general relativistic effects enter this disper- 
sion relation like the transverse electromagnetic waves 
by the ratios of uj/a. If the lapse function is one, 
the effect of gravity vanishes and correspond to the 
special relativistic version. 

9.1 Longitudinal Waves with — 3> kuQ S 
For — 3> kuo s , the dispersion relation given in Eq. 
(|51[) can be approximated as 
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The low frequency mode corresponding to ion acous- 
tic wave in the electron-ion plasma does not exist 
because the charge separation never occurs in that 



situation. For electron-positron plasma, we can ob- 
tain the dispersion law of the waves under consider- 
ation from Eq. ([52]) as 
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where X]j s = ^kBT os /4ire 2 no s is the the debye ra- 
dius. Equation ([59]) satisfies ^ 3> kuo s , if the second 
term in the right-hand side of this equation is small 
compared to the first term. This means that 



\ 2 Ds k 2 « 1, 



(54) 



which corresponds to rather long wavelength os- 
cillations k 2 « 2uj 2 s at any arbitrary relativistic 
plasma temperature. In the limit of zero gravity 
(i.e., when a — >■ 1) the above equation becomes 
uj 2 = 2uj 2 ps {\ + ^-X 2 Ds k 2 ), which is the SK [4] re- 
sult with 7 g = 1. 

9.2 Longitudinal Waves with — ks kuo s 
In this approximation the dispersion relation given 
in Eq. (IBTT) reduces to 
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(55) 



It follows from Eq. ([55]) with < uo s < 1 that in 
this case the branch of a periodically damped and 
growing (i.e., unstable) oscillations exist either for 
electron-positron or electron-ion plasma. When the 
plasma is freefalling onto the horizon Eq. ([55]) be- 
comes undefine, that is to say the isotropic pressure 
will not hold. For the electron-positron plasma Eq. 
([55]) can be written as 
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(56) 



Equation ([59]) is valid for long wavelength oscil- 
lations and correesponding to the same dispersion 
curve given in Eq. (155]) . 

9.3 Longitudinal Waves with — <C kuo s 
For - <C ku 0s , we obtain from the Eq. ([5T]) . 
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(57) 



This equation shows that a singularity take place 
for each fluid species at the point for which the 
fluid's velocity equals the half of the thermal ve- 
locity, Uq s = ^w^ s , i.e., the transonic radius be- 
gins to play a significant role for the longitudinal 
waves. The position of transonic radius of each fluid 
mainly dependent on their limiting horizon temper- 
ature and determines the temperature at any given 
radius. We have from Eq. (|57[) for electron-positron 
plasma 
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(58) 



This equation is valid for the range of large wave 
numbers 

X 2 Ds k 2 » 1. (59) 

Therefore, the formulae ([5Tf and (|57|) can be qual- 
itatively matched at Xjj s k 2 ~ 1 and can be consid- 
ered one and the same dispersion curve as various 
regions. 

10. Concluding Remarks 

We have made an analytical study of transverse and 
longitudinal waves in an ideal relativistic two-fluid 
plasma, to advance some results obtained by differ- 
ent authors and to include corrections due to effects 
of general relativity. 

Considering the transverse electromagnetic os- 
cillations the dispersion curve are described by dif- 
ferent analytical formulas. It follows from our anal- 
ysis that one purely real Alfven and high frequency 
electromagnetic modes are exist in high and low fre- 
quency limits. This results are is in agreement with 
the results of SK [10] for ultrarelativistic limit and 
of Daniel and Tajima [10] for general relativity. For 
very low or negligible frequency, damped and grow- 
ing modes exist and the same conclusion also follows 
from the numerical calculations of Buzzi et al [25j . 

For longitudinal oscillations, transonic radius be- 
gins to play a significant role. There exists only one 
real high frequency dispersion curve in high and low 
frequency limits as was found by SK in ultrarela- 
tivistic limit. For a very low frequency, damped 
and growing modes exist like the transverse electro- 
magnetic waves. The presence of damped modes 
demonstrates that energy is being drained from the 
waves by the gravitational field and growing modes 
point out clearly that the gravitational field is, in 
fact, feeding energy into the waves. 

I am grateful to the anonymous referee for point- 
ing out some mistakes and ambiguities in the first 
version of this manuscript and for giving me some 
references to improve the structure of this paper. 
I am glad to acknowledge the Editor Dr. Ronald 
C. Davision for sending me a paper related to this 
work. 
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